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Resonances in the reflection probability amplitude r(E) can occur in energy ranges in which 
the reflection probability R(E) = |r(_E)| 2 is 1. They occur as the phase <j){E) defined by r(E) = 
t*(E)/t(E) = le i2 ^ E ' undergoes a rapid change of ty radians. During this transition the phase 
angle exhibits a Lorentzian profile in that d<j>(E)/dE ~ — Eo) 2 + (h-y/2) 2 ]. The energy Eq 

identifies the location of a quasi-bound state, 7 measures the lifetime of this state, and t(E) is a 
matrix element of the transfer operator. Methods for computing and measuring these resonances 
are proposed. 



PACS numbers: 03.65.-w,42.50.Xa,85.25.Dq 

I. INTRODUCTION 

Quantum mechanical transmission resonances, as ob- 
served for example in the Ramsauer effect, have been 
known for a long time []]. At a transmission resonance 
the reflection probability is small, and may even be zero 
if the potential satisfies specific properties. Conversely, if 
the transmission probability is small the reflection prob- 
ability will be large. 

There is a large class of potentials for which the trans- 
mission probability is zero and therefore the reflection 
probability is one in certain energy ranges. In one di- 
mension these potentials have unequal asymptotic values 
on the left and right, with Vl < E < Vr. Here Vl 
(Vr) is the asymptotic potential on the left (right) and 
E is the energy of a particle, incident from the left. One 
such potential is shown in Fig. 1. This is a washboard 
potential, commonly encountered with biased Josephson 
junctions 0- Particles incident from the left with energy 
-10 < E < 19. leF will be reflected with 100% proba- 
bility, so that R(E) = |r(.E)| 2 = 1.0. It is the purpose 
of this work to show that the reflection probability am- 
plitude r(E) undergoes resonances in the sense that the 
reflection phase shift <f>(E), defined by r(E) = le t2 ^ E \ 
undergoes rapid phase shifts when the incident energy 
sweeps through a narrow energy range. These resonances 
correspond to energies at which quasistable states occur 
in the potential. The sharp resonances can be described 
by a Lorentzian line shape whose peak locates the energy 
of the metastable state and whose width is determined 
by the lifetime of the state. 

II. HAMILTONIAN 

We choose as hamiltonian 

n = £ + v{i) =Si +v ° cos(7) + Vi1 (i) 

Washboard potentials of this type are typically found in 
biased Josephson junctions Q|. In such cases the coor- 
dinate 7 is the phase difference across the junction, P is 



Washboard Potential 

V(x) = V_l x +V_0 cos(x/L) 
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FIG. 1: Washboard potential: V(x) = Vocos(x/L) + ViX, 
with Vo = -10 eV, Vi = 1 eV, L = 1 Nm, < x < 9.7 
Nm, and Vl = -10 eV, V R = 19.9 eV. H = P 2 /2M + V(x), 
M = 100m e . This potential has two metastable states, at 
E = 0.32 eV and E = 7.87 eV. 

the canonically conjugate momentum, M = Cjfi 2 /Ae 2 is 
expressed in terms of the junction capacitance Cj and 
fundamental constants e and h, Vq is the Josephson en- 
ergy, and V1/V0 = I/I c , where I is the bias current and 
I c is the critical current. For simplicity, we will inter- 
pret this hamiltonian as describing a particle of mass 
AI = 100m e with coordinate x = 7L in the potential 
whose parameters are as described in the caption of Fig. 
□ 



III. REFLECTION PHASE 

The phase shift can be computed by constructing the 
transfer matrix for the potential. This relates the input 
and output amplitudes on the left with those on the right 
(the notation is as used in 
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The reflection amplitude is r(E) — t 2 i(E)/tu(E). For 
reflecting boundary conditions these two matrix elements 
are complex conjugates of each other: 



t 21 {E)=t ll {E)* 



t n (E) = a(E) + ib(E) (3) 



In the neighborhood of a metastable state the the matrix 
element t n (E) has the form e-' e [(E ~ E ) - i(fry/2)], 
where the rapidly varying phase tan -1 (h , y/2/(E — Eq)) 
rides on top of a slowly varying 6(E). The total phase is 
4>(E) = 0(E) + tan" 1 ((hj/2)/(E - E )). From this we 
immediately derive a Lorentzian line shape: 
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When 0(E) is slowly varying over the narrow range 
of a resonance, the term dQ/dE can be neglected and 
d<p(E)/dE shows a Lorentzian lineshape. 

The real and imaginary parts of the matrix element 
tn(E) have been computed for the potential shown in 
Fig. n They are plotted in the neighborhood of the first 
excited resonance at E = +7.87 eV in Fig. |2] The phase 
shift (j)(E) computed from the real and imaginary parts 
of tn(E) are shown in Fig. [21 

The real and imaginary parts of these matrix elements 
can be written in the form 



a(E) = a(E)(E - E r ) 
b(E) = (3(E)(E-E l ) 



(5) 



They have zero crossings at E r and Ef. the sharper the 
resonance, the closer the crossings. The slopes a(E) and 
/3(E) are related to cos(#) and sin(0). Generally the 
slopes are only weakly dependent on E. When this is 
the case 



\t 11 (E)\ 2 ^\t 21 (E)\ 2 ^a(E) 2 +b(E)^ 
(a 2 +(3 2 )[(E~E Q ) 2 + (hj/2) 2 } 
Its inverse has a Lorentzian shape 
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The center Eo and halfwidth fry/2 of the peak are 
weighted averages of the zero crossings of the real and 
imaginary parts of tn(E): 
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Figure 01 plots the Lorentzian l/\tn(E)\ 2 in the neigh- 
borhood of a sharp resonance for the washboard poten- 
tial. Also plotted is the inverse of this function, which 



should be a simple parabola if the peak does in fact have 
a Lorentzian profile. When the slowly varying part of 
4>(E) is removed, so that (j)(E) varies rapidly through tt 
radians as E passes the resonance, a plot of d<fi(E)/dE 
vs. E produces a Lorentzian line shape centered at the 
resonance. The plots of l/\tn(E)\ 2 and d(f>(E)/dE are 
indistinguishable. 

Matrix Element, t_{ 1 1 }(E), Real & Imaginary Pts 
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FIG. 2: Real (negative slope) and imaginary (positive slope) 
parts a(E) and b(E) of the matrix element tn(E) — a(E) + 
ib(E). 



Reflection Phase Angle / Pi 




7.85 7.9 
Energy 



FIG. 3: Phase shift 4>(E) /n computed from the real and imag 
inary parts of the matrix element tu(E). 



IV. PHYSICAL INTERPRETATION 

Rapid phase changes in transmission and/or reflec- 
tion amplitudes are associated with the existence of 
metastable states. They have also received a nice inter- 
pretation by Wigner as time delays caused by these res- 
onances 0. That is, as a particle enters a region where 
a metastable state can exist, it undergoes multiple inter- 
nal reflections before exiting this region. These internal 
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FIG. 4: Reflection resonance for the washboard potential 
shown in Fig. 1. Also shown is its inverse, which should 
be a simple parabola if the resonance is a Lorentzian. 

reflections are responsible for the prolonged time delay 
before transmission or reflection. The Wigner time delay 
is measured quantitatively by St = Tid<j)/dE. This delay 
time has been studied under scattering conditions [SJ , re- 
flecting boundary conditions [|J, and metastable escape 
conditions 0,0- 

The phase shift calculated in Sect. Ill shows 
Lorentzian peaks near resonances riding on top of a 
slowly varying background. The slowly varying back- 
ground is due to the slight time delay for the reflected 
particle due to the longer distance it must travel before 
reaching the classically forbidden region. The fast vari- 
ation is the excess time delay spent in the well near the 
energy of the metastable (resonant) state. 

V. DIFFERENTIAL SPECTROSCOPY 

It may not always be easy to extract phase informa- 
tion from reflection amplitudes. An alternative method 
for locating reflection resonances would therefore be use- 
ful. One such method is illustrated in Fig. This is a 
schematic for a Differential Reflection Resonance Spec- 
trometer. A beam of thermal electrons (E ~ eV) is 
incident on a beam splitter from arm of the interfer- 
ometer. The beam is split into arms 1 and 2. Each of 
these beams is reflected off identical washboard poten- 
tials. The two reflecting potentials are biased at voltages 
— V and —V — SV 0- The separation SV between the 
two biasing potentials is chosen small, of the order of 10 
profile halfwidths. The potential V is scanned through 
the resonance, and the intensity profile of the interfering 
beams in arm 3 is recorded. The split beams undergo 
rapid 27r phase shifts approximately SV apart. As these 
two phase shifts occur the output intensity undergoes a 
rapid symmetric variation as shown in Fig. The sym- 
metry of this pattern and the separation of the rapidly 
fluctuating features by SV is a clear indication that a 



reflection resonance exists approximately where the first 
feature undergoes its most rapid variation. This result 
is insensitive to amplitude and phase inequalities in the 
two returning beams due to different arm lengths and 
unequal splitting in the beam splitter. 




FIG. 5: Reflection resonance interferometer. A beam of ther- 
mal electrons E ~ is incident on a beamsplitter. The two 
resulting beams are reflected off identical washboard poten- 
tials that are biased at —V and —V — 5V, with SV of the 
order of 10 halfwidths of the resonant peak. The bias poten- 
tial V is scanned through the resonance and the intensity of 
the recombined beam in the output arm is monitored. 



VI. LORENTZ PROFILES FROM INTENSITIES 

Line shape information can be extracted from intensity 
information. The total amplitude of the signal in the 
measurement arm of the interferometer is 

A(V) = a^niV) + a 2 e la2 r 2 (V + SV) (9) 

Here r\(V) — e l2 ^ v ^ is the reflection amplitude for par- 
ticles reflected from the potential biased at V in arm 1, 
and similarly for r 2 (V + SV) = e l2 ^ v+5v ^ in arm 2. The 
returning beams are weighted with complex amplitudes 
ciie tai , with cii > 0. The amplitudes dj are typically un- 
equal due to nonperfect beam splitting, and the difference 
in the angles on is related to differences in arm lengths. 
The observed intensity is 

I(V) = \A(V)\ 2 =aj + al + 

(10) 

2a ia2 cos [(ai - a 2 ) + 2(0(V) - <j>(V + SV))} 

A typical intensity profile is shown in Fig. [fj] For this 
pattern SV — 0.2eV, a\ = 1,cki = and a 2 = 0.7, a — 
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0.2 x 2n. The intensity varies between /Max and J m i n , 
where cos [*] assumes values +1 and —1. Several features 
of this profile are worth mentioning. (1) The pattern is 
symmetric about the midpoint V + \8V . This comes 
about because the intensity is a function of the differ- 
ence 2{4>{V) — 4>(V + SV)), and these two functions are 
identical except for the energy shift. Outside the range of 
rapid variation this difference is zero, while on transiting 
the resonance region the difference changes from to 2n 
radians, then back down to 0. (2) The pattern has five 
critical points. The two critical points near the resonance 
at — V = —7.9 eV occur as 2cf>(V) changes through 2n ra- 
dians and the two critial points near —V — SV occur as 
2<j)(V + SV) changes through 2tt radians. The fifth criti- 
cal point, at the symmetry point V + \SV , is due to the 
overlap of the shoulders of the two Lorentzians centered 
near V = — 7.9 eV and V = —8.1 eV, as will be shown 
shortly. 

The derivative of the intensity is 
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At /Max and / m j n where cos [*] = ±1, sin [*] 
passes through zero. As a result 2a\a 2 sin [*] = 
±yj (/Max — /)(/ — /min)- Therefore, the intensity data 
can be processed according to 
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The right hand side is the difference of two Lorentzians. 
This difference vanishes at the midpoint V + \5V , where 
the falling shoulder of the Lorentian around the res- 
onance near —7.9 eV is equal to the negative of the 
rising shoulder of the displaced resonance near —8.1 
eV (fifth critical point of dI(V)/dV). The expres- 
sion on the left of Equ.(ll) is difficult to compute 
near the critical points where / = /Max or / m i n . In 
these regions the intensity should be approximated in 
second order, and the singular expression simplifies to 
V \d 2 I{V)/dV 2 1 y/ (/Max - /min)/2. The absolute value of 
the right hand side of Equ.(ll) is plotted in Fig. 0for 
the intensity data shown in Fig. The absolute value 
consists of two Lorentzian profiles, one surrounding the 
resonance near —7.9 eV, the other around the displaced 
resonance near —8.1 eV. These two peaks are displaced 
vertically from the curve dI(V)/dV, shown below it. 



VII. CONCLUSIONS 

For scattering potentials possessing bound states, the 
phase of the matrix element tn(E) increases by ir radi- 
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FIG. 6: Intensity in the output arm as a function of bias- 
ing voltage V in the Differential Reflection Resonance Spec- 
troscopy setup shown in Fig. 5 (darker curve). For this 



simulation we used SV = 0.2 eV, r% = e 



i2<t>(V) 



and T2 = 



ae ie e i24,{v+sv) with a = Q 7 and = o.2 x 2tt. The plots of 
4>{V) and (f>(V + 5V) are shown in lighter curves after being 
shifted down by an equal phase. 
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FIG. 7: Plot of dl/dV for the intensity shown in 
Fig. |S| Vertically offset: plot of the function 
dl/dV/ \J (/Max — /)(/ — /min) , showing recovery of two 
Lorentzian line profiles at the appropriate biasing potentials. 



ans each time E increases through a bound state value 
(real pole of the ^-matrix) and increases rapidly through 
7r radians each time E passes through a transmission res- 
onance. For such resonances the transmission line shape 
l/\tn{E)\ 2 is well approximated by a Lorentzian. For re- 
flecting potentials of the type discussed here (c.f., Fig. QJ, 
the reflection amplitude r(E) — tn(E)* /tu(E) increases 
rapidly through 2?t radians each time E passes through a 
reflection resonance. For such resonances the line shape 
1/ |£n (Z?) | 2 is well approximated by a Lorentzian. The 
Lorentzian can be computed from the real and imagi- 
nary parts of the appropriate transfer or scattering ma- 
trix element (c.f., Equ.(5) and Fig. 0J. It can also be 
computed as d<f)(E)/dE (c.f., Fig. EJ). We have described 
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how this information can be extracted from unique in- 
tensity signatures obtained from a Differential Reflection 
Resonance Interferometer (c.f., Equ. (11) and Figs. 
and0). Quasibound states in washboard potentials are 
currently located experimentally by spectroscopic meth- 



ods [2| . This requires at least one of the states to be oc- 
cupied. The method proposed here provides a completely 
different approach for locating such states. Further, the 
states need not be occupied. 
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